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Abstract 

In the tunnel limit, the current-phase relation of Josephson junc- 
tions can be expanded as /(</>) = I\ sin</> + I 2 sin 20. Standard BCS 
theory predicts that I\Rn ~ A/e and I2/I1 ~ D, where R N is the 
resistance of the junction in the normal state, A is the superconduct- 
ing gap, and D -C 1 is the junction transparency. In the cuprates, 
the experimental value of IiRn (I2/I1) is much smaller (larger) than 
the BCS prediction. We argue that both peculiarities of the cuprates 
can be explained by postulating quantum fluctuations of the pairing 
symmetry. 

1 Introduction 

The proximity to the metal-insulator transition is well known to lead to 
an anomalous normal state of the high T c superconductors. It is therefore 
interesting to ask whether, apart from the d-wave symmetry of pairing, also 
the superconducting state of the cuprates can be regarded as unconventional. 
Most studies attempt to answer this question by considering the properties 
of quasiparticle excitations. For instance, photoemission experiments seem 
to support the conventional alternative, since sharp spectral functions have 
been observed at low temperatures O. On the other hand, the effect of 
strong correlations on the condensate has remained largely unexplored so 
far. This is surprising, since large quantum phase fluctuations are to be 
expected in the cuprates as a combined effect of the suppression of charge 
fluctuations and of the uncertainty principle. In particular, such fluctuations 
have been suggested to stabilize the RVB state and the presumably related 
pseudogap phase of the cuprates [§]. In this paper we argue that by a 
detailed analysis of the Josephson effect, new insights into the nature of 
quantum fluctuations in the cuprates can be obtained. 

Josephson junctions involving the cuprate superconductors have been 
studied mainly because they enable phase-sensitive tests of pairing symme- 
try ||. Attention has been paid especially to two particular types of Joseph- 
son junctions: grain boundary [|j and intrinsic || Josephson junctions. The 
best studied type of grain boundary Josephson junctions involves junctions 
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built on c-axis oriented films, where the weak link forms at the boundary of 
grains which are rotated around the [001] axis with respect to each other. 
Idealized junctions of this type are characterized by a planar interface and 
two angles 9\ and #2 between the interface normal and the crystallographic 
directions in the grains forming the junction. The properties of the junc- 
tions depend dominantly on the misorientation angle 9 = 62 — 0\. It is 
well known Q that the transparency D of grain boundary junctions decays 
exponentially with increasing misorientation angle 9, D{9) oc exp(— 9/9q). 
Thus, for 9 > #0 ~ 5°, grain boundary junctions are in the tunnel limit 
and (for 9 not too close to 45°, see Section 3) their current-phase relation 
I{4>) = I\ sin 4> + 12 sin 2<fi + . . . can be well approximated by I((f)) = I\ sin cp, 
neglecting the higher-order harmonics. 

2 The Josephson product 

A useful quantity characterizing the superconducting electrodes forming the 
Josephson junction is the product of the first harmonic I\ with the junc- 
tion resistance in the normal state, Rn- According to standard theory (for 
homogeneous featureless barriers) , this so-called Josephson product is inde- 
pendent of the junction area and of the barrier transparency, thus giving 
an intrinsic information about the superconducting banks. It is well known 
that at interfaces between d-wave superconductors, anomalous bound An- 
dreev levels may form ||. At temperatures larger than the energy of such 
Andreev levels, BCS-like theory for rough interfaces between d-wave super- 
conductors predicts 

(IiRn)bcs = (vr/4)(A/e) cos 29, (1) 

where A is the maximal superconducting gap. The measured Josephson 
product of cuprate grain boundary junctions 0] can be well described by 
I\Rn = & 2 (IiRn)bcs with a ^-independent renormalization factor a 2 ~ 
10 _1 . In addition to the 9 > 9q data of [|J], this functional form can be 
tested also for 9 = (which case can be realized in break junctions) with 
the result that average Josephson products of such junctions || are fully 
consistent with the grain boundary data. Moreover, in it has been shown 
that A is not depressed in the junction region, thus explicitly demonstrating 
the breakdown of the BCS prediction for I±Rn in the cuprates. 

There exists no generally accepted explanation of the small renormaliza- 
tion factor a 2 . One of the reasons is that the microstructure of Josephson 
junctions is typically quite complicated. In fact, it is well known that small 
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angle grain boundaries can be modelled by a sequence of edge dislocations, 
while at larger misorientation angles the dislocation cores start to overlap 
and no universal picture applies to the structure of the grain boundary. For 
large-angle grain boundaries, Halbritter has proposed |J that the junction 
can be thought of as a nearly impenetrable barrier with randomly placed 
highly conductive channels across it. If due to strong Coulomb repulsion 
only the normal current (and no supercurrent) is supported by these chan- 
nels, the small value of I\Rn follows quite naturally. 

In this paper we shall argue that the smallness of the Josephson product 
does not follow from the particular properties of the barrier, but is rather 
an intrinsic property of the cuprates. Such a point of view has been first 
advocated in ]l~0fl . However, that paper did not consider alternative more 
conventional explanations. In order to support our point of view, let us 
begin by discussing the Josephson product for intrinsic Josephson junctions 
in the c-axis direction. Such junctions can be viewed as an analogue of aft- 
plane break junctions (since the misorientation angle vanishes for both), but 
are preferable because of simpler geometry of the interface. Moreover, zero 
energy surface bound states which may develop at aft-plane surfaces because 
of the d-wave symmetry of the pairing state || do not form in the c-axis 
direction, simplifying the analysis of intrinsic Josephson junctions. 

Standard BCS theory applied to the case of tunneling between two- 
dimensional superconductors fl]J predicts (for coherent c-axis tunneling) 
that the c-axis critical Josephson current density is j% = (2e/h)N(0)(z) ( .t 2 ,) , 
where (...) denotes an average along the (two dimensional) Fermi line, z k 
is the wavefunction renormalization, and t k is the matrix element for c-axis 
tunneling between neighboring CuC-2 planes. -/V(O) is the bare (unrenor- 
malized) density of states, iV(0) = § dk(^ir 2 hvk)~ l , where the integration 
runs along the two dimensional Fermi line and v k is the bare Fermi velocity 
at the Fermi surface point k. Note that although we have assumed that 
the self-energy is only frequency dependent, for this geometry z k enters the 
expression for j% (and also the normal-state resistance, see below). 

Let us also note that the use of ordinary perturbation theory in deriving 



an expression for j c has been criticized recently ]12|| . However, our formula 
yields the correct answer for t/~ < A. This can be shown either by an explicit 
solution of a 4 x 4 Bogoliubov problem for two coupled planes with phase 
differences and ir between the planes, or by considering solutions to the gap 
equation in an infinite layered system with a finite c-axis total momentum 
of the Cooper pairs (for a similar calculation, see e.g. [|l3f| ) . 

The conductance per square in the normal state is given by Gn = 
(2e 2 /h)N(0)(z k t 2 k F^ 1 ) N , where T k is the inverse lifetime of the quasipar- 
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tides and the index N in the Fermi line average means that the quan- 
tities are to be evaluated in the normal state. Therefore standard the- 
ory predicts for the Josephson product of intrinsic Josephson junctions 
I\Rn = jiG]^ 1 ~ e~ l (zkt\) / (zkt\T~^ 1 ) n . In conventional superconductors 
Rn can be measured at low temperatures in a sufficiently large magnetic 
field. This is impossible for the cuprates and thus Rn is usually defined as 
the c-axis resistivity at T c . 

Unfortunately, due to the unknown temperature dependence of z k , the 
theoretical I\Rn can not be directly tested by experiment. In order to 
overcome this problem, in M] instead of the usual Josephson product a re- 
lated characteristic of intrinsic Josephson junctions has been studied, namely 
the product of the critical current I\ and of the resistivity R$ in the re- 
sistive mode of the junction (at low temperatures). Since the conduc- 
tance per square in the resistive mode of the superconductor is [14] Gs ~ 
(8e 2 /h)N(0)(z k tl) node /A, standard BCS-like theory predicts hR s = jiGg 1 ~ 
(tt /2)(A/ e)(z k tl) / (z k tl) node . In Q, hR s ~ A/e has been found experi- 
mentally and good agreement with theory has been claimed, since momentum- 
independent t k and Zk were assumed. However, according to band structure 
calculations |L5|], t k is strongly suppressed in the nodal directions. If this 
modulation of the tunnel matrix element t k is taken into account and the 
presumably only moderate A;-space dependence of z k is neglected, the exper- 
imental I±Rs is seen to be drastically reduced with respect to the theoretical 
predictions. 

Thus we have shown that although the barriers in grain boundaries and 
in intrinsic Josephson junctions are of very different nature, both types of 
junctions exhibit a suppressed Josephson product. Therefore we believe that 
this suppression is not due to specific barrier properties as suggested in || , 
but rather due to some intrinsic property of the high-T c superconductors. 



3 The second harmonic of the current-phase rela- 
tion 

Since the second harmonic I2 is not forced by symmetry to depend on the 
angles its Josephson product can be estimated (for temperatures larger 
than the energy of anomalous Andreev levels) using the standard BCS theory 
as I2RN ~ DA/e. Comparison with Eq. (|l]) implies that for junctions with 
9 ~ 45° the d-wave symmetry of pairing leads to a suppression of 1%, and 
I2 may become comparable to I\. This has in fact been observed in two 
different types of 45° grain boundary Josephson junctions [16, 17]. 



4 



However, the results of [16, 17] are quite mysterious, if we take into 



account the actual experimental setup. In fact, standard BCS theory with 
ideal featureless barriers implies that in order that I-Z2I ~ \Ii\, the average 
misorientation angle 9 would have to be given with a precision ~ D, where 
D ~ 10~ 3 . This is not realistic and therefore two alternative explanations 
have been proposed, in both of which the origin of the anomalously large 
l^/^il ratio has been sought in the barrier properties. 

(i) Faceted scenario has been considered as an alternative explanation 
for symmetric 45° junctions (i.e. junctions with nominal geometry 6\ = 0° 
and 62 = 45°), in which I/2I > I h I has been found [16|. It takes into 
account the faceting of the grain boundary and also the twinned nature of 
the (orthorhombic) YBCO thin films. Due to both of these features, the 
junction can be viewed as a parallel set of and ir junctions It has 

been shown [19| that in such a case spontaneous currents are generated 
along the interface, the ground state energy of the junction is minimized at 
a macroscopic phase difference ±7r/2, and consequently the current-phase 
relation is dominated by the second harmonic Ii- 

In what follows we analyze quantitatively whether the faceted scenario 
can apply to the results of [16|. Let us denote the current densities cor- 
responding to the harmonics Ii (with i = 1,2) as ji and introduce the 
Josephson penetration depth of the junction, Aj = ($o/47rA//oj2) 1//2 - More- 
over, let the local critical current density in the and ir junctions be ±jo, 
their typical length a, and the bulk penetration depth be A. Then, since 
a ~ 0.01 — 0.1/um, A ~ 0.15/iin, and Aj ~ 3fim (estimated making use of 
[@] J2 ~ 10 4 A/cm 2 ), the inequalities 7rA 3> a and A 2 3> a A are well sat- 
isfied. Following [|H]] it is easy to show that these inequalities guarantee 
that the spontaneously generated currents along the Josephson junction can 
be calculated within perturbation theory. If we denote the total junction 
length by L, then a straightforward calculation yields J2/J0 ~ MoJoaA 2 /<I>o 
and ji/jo ~ ( a /L) 1 ^ 2 . The equation for j\ is a random walk-type formula, 
indicating that j% averages to zero in a sufficiently long junction. After 



some algebra the above equations are seen to imply j'2/ii ~ yLA/(47rAj). 
Therefore, standard theory predicts that I2 > Ii can be realized only in 
sufficiently long junctions with L > 47rAj/A. This requires L > 500/mi, 
whereas in ]l6| much shorter junctions with L ~ 1/mi were studied. 

(ii) Pinhole scenario has been proposed in [17] as an explanation for 
symmetric 45° junctions (i.e. junctions nominally characterized by 82 = 
—9\ = 22.5°). It views the barrier as basically impenetrable, the conduction 
being due to randomly placed highly conductive pinholes. This explains 
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quite naturally the small value of the effective barrier transmission and, at 
the same time, the large value of \I2/I\\- Note that in order to explain the 
small value of the Josephson product, in addition to pinholes also Halbrit- 
ter's conductive channels || have to be postulated, which are assumed to 
be highly conductive only in the normal and not in the superconducting 
channel. Because of this ad hoc nature of the pinhole scenario, and mainly 
because of the absence of higher harmonics in /(</>) at 4 K jl7| whose presence 
it predicts, we believe that the pinhole picture should be discarded. 

Thus we conclude that the large value of the second harmonic I2 (com- 
pared with predictions of the standard BCS theory) is most probably not 
an extrinsic (barrier-related) effect, but rather an intrinsic property of the 
cuprates. 



4 Microscopic implications 

The two apparently unrelated experimental facts, namely the suppressed 
Josephson product I±Rn and the enhanced ratio \l2/h\, can be explained 
by a single assumption that in the cuprates some mechanism is operative 
which leads to a suppression of I\, while leaving Rn and I2 intact. In what 
follows we describe one such mechanism which we believe to be the most 
promising one. Namely, we suggest that at low temperatures the super- 
conducting state of the cuprates supports fluctuations of pairing symmetry 
towards s-wave pairing (which pairing is expected to be locally stable within 
several microscopic models of the cuprates). Such fluctuations presumably 
do not affect Rn, while they do influence the Josephson current. In simplest 
terms, if we denote the phases of the superconducting grains forming the 
junction as fa, then the fluctuations renormalize the first and second har- 
monics by the factors (e 1 ^ 1 ) (e 1 ^ 2 ) and (e 2l<?il )(e 21 ^ 2 ), respectively, where (. . .) 
denotes a ground-state expectation value. Thus experiment requires that the 
fluctuations have to be of such type that |(e*^)| = a ~ 0.3 and |(e 2j<?!> )| ~ 1. 
Precisely this behavior is expected if the d-wave order parameter fluctuates 
towards s-wave pairing. 

Now we proceed by introducing a minimal model of such fluctuations. 



Unlike in standard literature on this subject (see, e.g., [20| and references 
therein), we assign a phase field (fi to each bond i of the square Cu lattice. 
In other words, cpi lives on the sites of the dual lattice. Since large phase 
fluctuations are expected, the compactness of phase fluctuations is explicitly 
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Figure 1: Ground state of the Heisenberg model on an elementary plaquette. 



taken into account and the model reads 

H = flfluct + [-Vcos(2pi - 2tpj) + Wcos&i - ipj)] , (2) 

where Ha uct = —m~ 1 J2id 2 /d l Pi- The second term in Eq. @ (with 
denoting a pair of nearest neighbor sites) describes for < W < V a super- 
conductor with a dominant d-wave pairing V + W and subdominant s-wave 
pairing V — W. Since we concentrate on the q ~ (ir, it) fluctuations from d 
to s-wave pairing, as a first approximation we do not take into account the 
coupling of phase fluctuations to electromagnetism. 

In order to gain insight into the microscopies of H^ uct , let us recall that 
the RVB fluctuations in an elementary plaquette (in a hole-free region) favor 
a minus sign between the two different valence bond configurations (see 
Fig. 1). If these configurations are thought of in terms of Bose condensates 
of valence bonds, this means that energy is gained for a relative phase of 
the x and y condensates arg(^~ 1) = ±vr/2. RVB processes are thus seen 
to frustrate the phase ordering dictated by the second term in Eq. (|2|). 

Since the other singlet state of the elementary plaquette (with a plus sign 
between the valence bond configurations) lies higher in energy, the effective 
Hamiltonian in the singlet sector is -£Trvb °c A] Al +x+y Ai +x Ai +y + H.C., 
where A] creates a singlet on bond i. Neglecting the amplitude fluctuations 
of Aj, we can write H R vb = JJ2i cos(^ - cp i+x + (pi +x+y - <Pi+ v ). In what 
follows we replace -fffl uc t by -ffRVB in the quantum model Eq. (|2|), thereby 
obtaining a completely classical toy model. 

The toy model permits a simple mean field analysis: We consider two 
types of solution, both of which live on two sublattices A and B, and 
we explicitly disregard four-sublattice solutions, since they correspond to 
translation-symmetry breaking states on the original Cu lattice. In the first 
type of solution ipA = <p and (fs = and there are two possibilities: either 
(p = 7T (cZ-wave), or ip < it (a complex mixture of d and s, to be called d+is). 
In the second type of solution we assume a disordered state of such type that 
on sublattice A, the phase equals and ir with probabilities (1 + a)/ 2 and 
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Figure 2: Mean field phase diagram of the toy model. 



(1 — oi)/2, respectively, and on sublattice B the values and tt are inter- 
changed. In this case the macroscopic symmetry is of the d-wave type, with 
renormalized averages (e i(pA ) = -(e i<fiB ) = a and (e 2itpA ) = (e 2iipB ) = 1. 
Minimization of energy (with respect to ip or a) leads to the phase diagram 
shown in Fig. 2. The renormalization factor in the disordered d-wave state 
(which presumably corresponds to a homogeneous but strongly fluctuating 
d-wave phase in the model Eq. (^)) is a = (W/J) 1 / 2 . Thus the cuprates 
correspond to the region W ~ J/10 and J < AV in Fig. 2. 



5 Conclusions 



The anomalous effects observed in cuprate grain boundary and intrinsic 
Josephson junctions can be explained by a single assumption of strong quan- 
tum phase fluctuations at q = (vr,7r), which leads to |(e^)| = a ~ 0.3 and 
|(e 2l< ^)| ~ 1. In addition, our picture implies that the Josephson product for 
junctions between the cuprates and low-T c superconductors is renormalized 
by the factor a, in semiquantitative agreement with experiment [21]. It also 
may be relevant for the experiment [22], where a large second harmonic has 
been found in a c-axis Josephson junction between YBCO and Nb. 
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